PROCEEDINGS 

°F SCIENCE 



Charm semileptonic decays and \Vcs(d)\ from heavy 
clover quarks and 2+1 flavor asqtad staggered 
ensembles 



Jon A. BaileyV D. Du,^ A.X. El-Khadra,^ Steven Gottlieb,^ 
A.S. Kronfeld,^ R.S. Van de Water,^ and R. Zhou' 

Email: 



R.D. Jain, 



jabsnuggmail . com 



"Department of Physics and Astronomy, Seoul National University, Seoul, 151-747, ROK 
^Physics Department, University of Illinois, Urbana, IL 61801, USA 
'^Department of Physics, Indiana University, Bloomington, IN 47405, USA 
'^Theoretical Physics Department, Fermilab, Batavia, IL 60510, USA 

Fermilab Lattice and M\LC Coiiaborations 

By combining experimentally measured partial branching fractions for the semileptonic decays 
D — > K£v and D n£v with lattice calculations of the form factors f'^^^iq^) and f'^'^^iq^), one 
can extract the CKM matrix elements \Vcs\ and \Vcd\- We are calculating the form factors by using 
Fermilab charm and asqtad staggered light and strange quarks on 2+1 flavor asqtad staggered 
ensembles generated by the MILC Collaboration. We vary the light valence quark masses from 
0.4ms to 0.05ms (m, is the strange sea-quark mass), and the lattice spacings, from about 0.12 fm 
to about 0.045 fm. We extrapolate to the physical light-quark mass and the continuum limit 
using heavy-light meson staggered chiral perturbation theory about the SU(2) and SU(3) limits, 
compare the resulting (preliminary) form factor shapes with experiment, and discuss our errors. 



The 30th International Symposium on Lattice Field Theory 
June 24 - 29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



D K(7t)£v and \Vcs{d) I from heavy clover on 2+1 flavor asqtad MILC ensembles 



Jon A. Bailey 



1. Introduction 

Tests of second row and column unitarity are limited by the uncertainty in \Vcs\, determined 
from (semi)leptonic decays, while assuming three-generation unitarity leads to the very precise 
(Standard Model) values for |Vci(rf)| [|l|]. Direct determinations of |V„(rf)| are being improved by 
combining increasingly precise lattice calculations of decay constants and form factors ^] and 
information from experiment [Q, |5|]. The results have been used to improve tests of unitarity and, 
by comparing with the Standard Model values, to validate lattice calculations of B decay constants 
and form factors. Ab initio calculations of hadronic weak matrix elements of B decays are central 
in a number of ongoing searches for new physics [^. 

In the rest frame of the D meson, the semileptonic branching fractions and CKM matrix ele- 
ments are related by 

d^^ 24^'^^^-' 'P^' ' ^^-^^ 

where we neglect the lepton masses, x = s,d denotes the strange or down quark, P = K,n is the 
daughter meson, = {pd—PpY is the invariant mass of the leptons, and the form factor ff^^{q^) 
is defined in terms of the hadronic matrix element of the flavor-changing vector current ^ = ixy^c: 

{p\r\D) = U{q^) [PD+PP- y +Mq ) ^2 ^ - (1-2) 

Given the partial branching fractions, dT/dq^, from experiment and the form factor normalizations. 



e-g; f+^^{q^ = 0), from theory, one can use Eq. ( pTll) to extract the CKM matrix elements \Vcx\- 

Below we briefly describe our method and data set before focusing on the current status of the 
chiral-continuum extrapolation, preliminary results for the form factor shapes, and our anticipated 
errors. At present the analyses are blinded; that is, the absolute normalizations of the form factors 
are hidden. 



2. Method 



For lattice calculations a convenient parametrization of the form factors is 

{P\V\D) = .Jl^^{vf\\{Ep)+p^f^{Ep)) , (2.1) 

where v = po/tno is the 4- velocity of the D meson, p± = pp — {pp ■ v)v, Ep = (m^ + mp — 
q^)/{2mj}) is the energy of the recoiling P meson, and is the lattice version of the flavor- 
changing vector current 

We construct the lattice current out of light staggered and heavy clover fields |^. To match 
the continuum currents, the lattice current must be properly normalized. The matching factors Z^n 
nearly equal the geometric mean of degenerate vector current renormalization factors Zy4 and Zy* , 
which we are calculating nonperturbatively [^. The correction factors Py/j are constructed to be 
close to unity, and a subset of our collaboration has calculated them in one-loop lattice perturbation 
theory if]: 

{P\r^\D)=Zy,{P\V,^,\D), Zy,=py,^Zy^Zy,, p^, = 1 + ^ ( ) • (2-2) 
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Although the deviation of PyM from one is formally ^{ag), its numerical coefficient is small, and 
the ^(tts) correction, very small. To protect against analyst bias, we blind the analyses by intro- 
ducing multiplicative offsets in Pyu. 

The matrix elements (P|Vct|D) can be obtained from 3-point and 2-point correlators of the P 
and D mesons. We extract {P\'V^x\D) from the correlator ratio 

1 cgM^(f,r;pp) Ep I le-'"'-' 

where ^p^ = (1, pp), and C2 are 2-point P- and D-meson correlators, respectively, and Cf^^ 
is a 3-point correlator in which the flavor-changing vector current destroys the D meson at rest 
and creates the P meson with momentum pp. T and t are the source-sink separation and current 
insertion time in C3, respectively, and C3, C2 are averages of 3-point and 2-point correlators; the 
averages were designed to suppress oscillations from opposite-parity states [^]. For insertion times 
t far from source and sink, \<^t <^T, the temporal and spatial components of the ratio respectively 
approach the form factors /n and /l. 



3. Data 



We analyze approximately unitary data generated on MILC asqtad ensembles with several sea- 
quark masses > O-OSm^ and four lattice spacings < 0.12 fm (Table [l]). The f -meson 2-point cor- 
relators and 3-point correlators are generated at momenta pp = (0,0,0), (1,0,0), (1, 1,0), (1,1,1), 
(2,0,0), in units of 2n/aL. The D-meson 2-point correlators are generated with the D meson at 
rest. We use local interpolating operators for the P 2-points, and smear the D interpolating oper- 
ators in the 2-point and 3-point correlators with a IS charmonium wavef unction [Q]. To increase 
statistics, the correlators are generated at four source times spaced evenly along the time direction. 
Autocorrelations are suppressed by randomizing the spatial location of the source with configura- 
tion. We generate the 3-point correlators at four source-sink separations T to construct the averages 
C3 (Eqs. (37,38) of Ref. [^) and to minimize errors while avoiding excited state contamination [jlcj]. 



4. Chiral-continuum extrapolation 

For sufficiently small energies, quark masses, and lattice spacings, S;^PT describes the energy, 
quark-mass, and lattice-spacing dependence of the form factor data in a model-independent way. 
Below we show fits of data through pp = (1, 1,0) to SU(3) and SU(2) S;tPT. The key to the plots is 
in Table [l[ The fits in Figs. |l]-0 are fits to the data for /l and f\\ separately; we have not imposed 
the kinematic constraint at = 0. Restricting the data in the fit to have momenta less than (1,1,1) 
ensures reasonable behavior in the chiral expansion even on the coarse OArUs ensemble, where the 
corresponding recoil energies are greatest. 



We obtain the SU(2) S^PT fit functions by integrating out the strange quark in SU(3) S^PT 1 1 1 



12]. We include the strange-quark mass dependence in the SU(2) S;^PT LECs by expanding in the 
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Figure 1: Fits of /^^'^ data to SU(3) (left) and SU(2) (right) S^PT. The fit functions include the chiral 
logarithms and analytic terms at NLO and analytic terms at NNLO. Errors are statistical, from bootstrap 
ensembles. 





r\EK nEx 



Figure 3: Fits of /f data to SU(3) (left) and SU(2) (right) S;t;PT. The fit functions include the chiral 
logarithms and analytic terms at NLO and analytic terms at NNLO. Errors are statistical, from bootstrap 
ensembles. 



4 



D K(7t)£v and \Vcs{d) I from heavy clover on 2+1 flavor asqtad MILC ensembles 



Jon A. Bailey 



a (fm) 


X n, 


amf'/amf'' 


Key 


A'^conf 






0.12 


20^x64 


0.02/0.05 


n 


2052 


0.02/0.0349 


0.1259 




20^x64 


0.01/0.05 


□ 


2259 


0.01/0.0349 


0.1254 




20^x64 


0.007/0.05 


□ 


2110 


0.007/0.0349 


0.1254 




24^x64 


0.005/0.05 


□ 


2099 


0.005/0.0349 


0.1254 


0.09 


28^x96 


0.0124/0.031 


o 


1996 


0.0124/0.0261 


0.1277 




28^x96 


0.0062/0.031 


o 


1931 


0.0062/0.0261 


0.1276 




32^x96 


0.00465/0.031 


o 


984 


0.00465/0.0261 


0.1275 




40^x96 


0.0031/0.031 


o 


1015 


0.0031/0.0261 


0.1275 




64^x96 


0.00155/0.031 


o 


791 


0.00155/0.0261 


0.1275 


0.06 


48^x144 


0.0072/0.018 


A 


593 


0.0072/0.0188 


0.1295 




48^x144 


0.0036/0.018 


A 


673 


0.0036/0.0188 


0.1296 




56^x144 


0.0025/0.018 


A 


801 


0.0025/0.0188 


0.1296 




64^ X 144 


0.0018/0.018 


A 


827 


0.0018/0.0188 


0.1296 


0.045 


64^ X 192 


0.0028/0.014 


V 


801 


0.0028/0.0130 


0.1310 



Table 1: Ensembles, light valence masses, and hopping parameters. The up-down valence mass is equal to 
the up-down sea mass, while the strange valence mass is approximately physical, tuned via the K mass. The 
charm mass is also approximately physical, tuned via the Ds mass. "Key" is the legend for plots in Sec. 4. 




Figure 4: Fits of /jj^^^ data to SU(3) (left) and SU(2) (right) S^PT. The fit functions include the chiral 
logarithms and analytic terms at NLO and analytic terms at NNLO. Errors are statistical, from bootstrap 
ensembles. 

deviation from the physical strange-quark mass. The NLO terms absorb the scale dependence of 
the chiral logarithms, and we also include the analytic NNLO terms. 

The SU(2) fits above have greater p-values than the SU(3) fits; the difference is pronounced 
for f^^^ (Fig. ^. However, the fit results are not obviously superior, and in some cases the chiral- 
continuum extrapolated values differ by a few percent. 

The fits to f^^^ have somewhat small p-values. At higher energies, the errors of the chiral- 
continuum extrapolated curves for ff^^"^ grow, and the curves for f^^'^ and f^^'^ exhibit appar- 
ently unphysical inflection. These behaviors at higher energies may reflect the absence of data 
from the more chiral and finer ensembles; adding data at higher momenta from these ensembles 
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may decrease the errors and eliminate this inflection. 

The factors Zy4 and Zyj are preliminary, and we are updating them. Before the final fits, the 
data must be shifted to the retuned Kc values. We have generated additional data on one of the 
coarse ensembles to correct for the error in jq. -tuning and estimate the remaining systematic error 
due to uncertainty in the value of the (retuned) Kc- 



5. Results 

Due to suppression by the heavy quark mass, the form factors f^^^ are dominated by the 
form factors /f for which our fits are very well-behaved. By normalizing the form factors to 
convenient fiducial points, we compare the shapes obtained from LQCD and experiments; this ap- 
proach eliminates the need for any assumption about the normalization of the experimental data. 
Below we overlay fiducially normalized ;^PT curves from our fits (S;;^PT extrapolated to the phys- 
ical light quark mass and continuum limit) and form factor shapes from CLEO and BABAR [Q, ^. 
The errors on the experimental (blue and violet) data points are from the full covariance matrix. 




Figure 5: Preliminary LQCD results for the shape of f+^"{q^), with statistical errors only, compare fa- 
vorably with CLEO-c data, with total errors . On the left the curve is the chiral-continuum extrapolated 
shape from SU(3) S;|fPT; on the right the curve is from SU(2) S^PT. The curves agree within statistics. 



including systematics. Even though we omit systematic uncertainties in the lattice results in the 
above plots, the qualitative agreement between the curves and experiment is perfectly acceptable 
for both SU(3) and SU(2) ;^PT. Quantitative tests can be performed by fitting the lattice results and 
experimental data separately to the z-expansion [13]. Once the quantitative compatibility of the 
lattice and experiment form factor shapes is verified, simultaneously fitting the lattice results and 
experimental data will yield the CKM matrix elements IV^ I and \Vcd\- 

From the SU(2) (SU(3)) fits above, the statistical errors in /+^'^(0) are 4% (4%), and those in 
/+^^(0)> 2.3% (3%). Important systematic errors are from heavy-quark lattice artifacts, the error 
in ri, and the error in the axial couphng g,^. Naively updating the gji error reduces the projected 
systematics to 3.4% [0, 15 1. A careful estimate of all systematics reflecting the entire data set 
has yet to be made; the difference between our present SU(2) and SU(3) curves is in some cases 
comparable to the other errors. 

Fermilab is operated under contract DE-AC02-07CH11359 with the U.S. DOE; J.A.B. is 
supported by the Creative Research Initiatives Program (2012-0000241) of the NRF grant funded 
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Figure 6: Preliminary LQCD results for the shape of f+^'^iq^), with statistical errors only, compare favor- 
ably with CLEO-c [^] and BABAR data with total errors. On the left the curve is the chiral-continuum 
extrapolated shape from SU(3) S^PT; on the right the curve is from SU(2) S;)fPT. The curves agree within 
statistics. 



by the Korean government (MEST); this work was supported by grants DE-FG02-91ER40661 
(S.G., R.Z.) and DE-FG02-91ER40677 (D.D., R.D.J., A.X.K.) of the U.S. DOE. 
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